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Abstract 

We give a non-perturbative proof of a gradient formula for beta functions of two-dimensional quantum 
field theories. The gradient formula has the form diC — —{gij + ^9ij + bij)l3^ where /?-' are the beta func- 
tions, c and Qij are the Zamolodchikov c-function and metric, bij is an antisymmetric tensor introduced 
by H. Osborn and Agij is a certain metric correction. The formula is derived under the assumption 
of stress-energy conservation and certain conditions on the infrared behaviour the most significant of 
which is the condition that the large distance limit of the field theory does not exhibit spontaneously 
broken global conformal symmetry. Being specialized to non-linear sigma models this formula implies a 
one-to-one correspondence between renormalization group fixed points and critical points of c. 



Contents 

1 Introduction 

2 The general gradient formula 

3 Zamolodchikov's formula 

4 The proof of the gradient formula (first steps) 

5 Sources and operations 

5.1 General formalism 

5.2 Summary 

6 The Callan-Symanzik equations 

7 The proof continued 

7.1 The IR condition and the sum-rule 

7.2 The term rfg 

7.3 The identity d^J^^x) = p^d^J^{x) 

7.4 E4 is an IR error term 

8 Conclusion of the proof 

9 Discussion 

9.1 Contact term ambiguities and scale dependence 

9.2 The infrared condition: an example 

9.3 Bare gradient formula 

9.4 Dressing transformations 

9.5 Renormalization group transformation as a flow of couplings 

10 Final comments 



1 



1 Introduction 



Change of scale in quantum field theories (QFTs) is governed by renormahzation group (RG) 
transformations. If a space of theories is parameterized by coupling constants {A*} the RG 
transformations are governed by a beta-function vector field 



for some metric G^^{X) and potential function S{X) defined on the theory space, has some 
history. One of the earliest papers devoted to this question was [2j. It was suggested in that 
paper that RG flows are gradient flows in a wide variety of situations. Gradient flows have 
some special properties. Thus, if the metric G^^ is positive definite, the scale derivative of 
the potential function is negative definite 



and therefore S monotonically decreases along the flow. This demonstrates irreversibility of 
the RG flows and forbids limiting cycle behaviour. Another appealing property of gradient 
flows is that the matrix of anomalous dimensions dijS^ is symmetric and thus their eigenvalues 
at critical points, that give critical exponents, are always real. 

The flrst perturbative computations in support of this idea were done for four- dimensional 
theories [3] . Later more evidence was found in the context of two dimensional general sigma 
models [1] , [5] . In [5] a gradient formula of the form (11. 2p was formulated for such models and 
shown to hold up to two loops for a particular class of sigma models. A crucial ingredient for 
a gradient formula for general sigma models was the introduction of the dilaton fleld [6], [7j. 
It was shown in [8] , [10] that including the dilaton couplings into a general sigma model one 
flnds that the vanishing beta function equations are equivalent to critical points of a certain 
functional at the leading order in a'. A gradient formula of the form (11.21) was checked for 
general sigma models in [9] to the flrst two orders in a'. In string theory conformal sigma 
models describe strings propagating on the sigma model target manifolds. The sigma model 
couplings parameterize a metric Gjj, antisymmetric tensor Bjj and a dilaton fleld $ deflned 
on the target space manifold. The gradient property (II. 2p attains a special signiflcance in this 
context becoming a manifestation of the string action principle. The condition for conformal 
invariance is that the beta functions vanish: (3^ = (3^ = /5* = 0. It is equivalent to string 
equations of motion. The gradient property (II. 2p thus means that the string equations of 
motion arise by varying a functional of couplings - which can be identifled with the string 
action functional. 

Another reinforcement of the gradient conjecture (11.21) for two-dimensional theories came 
from the Zamolodchikov c-theorem [TT]. The last one is a general theorem applicable to 
unitary 2D theories that states that there is a function c on the space of theories that 
monotonically decreases along the RG flows and coincides with the Virasoro central charge 
at flxed points. (We give a slightly modifled proof of this theorem in section [3]). The theorem 




= nx) . 

•adient flows, that is 



(1.1) 



(1.2) 



(1.3) 
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was proved by constructing c whose scale derivative takes the form of the right hand side 
of fll.3p with a certain positive definite metric. It was natural to conjecture that a gradient 
formula of the form (11.21) holds with S being the c function and Gij - the Zamolodchikov 
metric. This was shown to hold at the leading order in conformal perturbation theory near 
fixed points [11], [13]. In the context of nonlinear sigma models this idea was discussed in |15j . 
It was argued in [15] that for the purposes of string theory the c-function cannot provide 
a suitable potential function (we comment more on this in section [TOl) . Other potential 
functions for RG flows of nonlinear sigma models were considered in [i2\, ^15], [H], [T7], [TS] 
which were shown to be related to the central charge and to each other. In [IS] a potential 
function for nonlinear sigma models was constructed assuming the existence of a sigma model 
zero mode integration measure with certain properties. It was shown that a measure with 
the required properties can be constructed infinitesimally but a proof of the integrability 
of that construction is still lacking. An essential tool proposed in [18] for deriving gradient 
formulas was the use of Wess-Zumino consistency conditions on local Weyl transformations 
in the presence of curved metric and sources. This technique was applied in [19] to a class 
of quantum field theories subject to certain power counting restrictions. It was shown that 
for these theories a gradient formula holds of a slightly different form than (11.21) : 

^^c = -g,,(3' - h,,(3^ (1.4) 

where c and g^j are the Zamolodchikov's metric and c-function [TT] and 6^- is a certain 
antisymmetric tensor. The necessity to introduce an antisymmetric tensor along with the 
Zamolodchikov's metric can be demonstrated by the use of conformal perturbation theory. 
Thus it was shown in [23] by explicit perturbative calculations that the one-form Qijd^c is 
not closed for some fiowqj. Still, as we will explain in the next section, Osborn's gradient 
formula (II. 4p . although very inspiring, falls short of providing a general gradient formula. 
The main content of the present work is a derivation of a gradient formula that generalizes 
formula (ll.4p to a much wider class of theories that includes nonlinear sigma models as well. 

To finish the historical overview we mention here that a general gradient formula was 
proven for boundary renormalization group flows in two dimensions [22] ■ Such flows happen 
in QFTs defined on a half plane (or a cylinder) when the bulk theory is conformal but 
the boundary condition breaks the conformal invariance. One of the implications of the 
boundary gradient formula is a proof of Affleck and Ludwig's (7-theorem ^21j which is a 
statement analogous to Zamolodchikov's c-theorem. A string theory interpretation of this 
gradient formula is that it provides an off-shell action for open strings. The boundary 
gradient formula was proved under certain assumptions on the UV behaviour which are 
reminiscent of the power counting restrictions of [19j. Nevertheless we will show in the 
present paper that any assumptions of this kind can be dispensed with in proving a bulk 
gradient formula. 

The paper is organized as follows. In section 2 after introducing some notations we 
explain in more detail Osborn's gradient formula (ll.4p and the assumptions that went into 
proving it. We then state our main result - a general gradient formula (I2.13P and discuss 
the assumptions needed to prove it. In section 3 we give a proof of Zamolodchikov's formula 
and recast it in the form that we use as a starting point for proving the gradient formula. In 
section 4 the first steps of the proof are explained. At the end of those steps we express the 
quantity diC + gijP^ + bijP^ built from the elements present in (ll.4p via three point functions 

^The obstruction to closedness occurs at the next to leading order in perturbation. 
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with a certain contact operator present in them. To analyze these three point functions we 
develop a sources and operations formalism in section 5. A short summary of the formalism 
is provided in subsection 5.2. After discussing the Callan-Symanzik equations in section 6 
we resume the proof in section 7 putting to use the Wess-Zumino consistency conditions 
on the local renormalization operation and our infrared assumptions. At the end of section 
7 an infrared regulated gradient formula is obtained. In section 8 the proof is concluded 
by removing the infrared cutoff. Section 9 contains a discussion of the properties of the 
gradient formula and the assumptions used in proving it. In section 9.5 the gradient formula 
is specialized to the nonlinear sigma model case and a proof is given of the correspondence 
between RG fixed points and stationary points of c. In section 10 we conclude with some 
final remarks. 



2 The general gradient formula 

In this paper we consider two-dimensional Euclidean quantum field theories equipped with 
a conserved stress-energy tensor T^y{x). The stress-energy tensor measures the response of 
the theory to metric perturbations, so that if Z[g^y\ is a partition function defined on a 
2-dimensional plane with metric g^u{x) = + 5g^y 

5\tiZ = ]- [[d^x{6g^,T'^''{x)). (2.1) 



2 . _ 

In two dimensions any metric can be made conformally flat so that g^u{x) = ^'^{x)5^y where 
the function sets the local scale. A change of local scale is generated by the trace of 
stress-energy tensor 6(x) = g^^T^i,{x) 

For correlation functions computed on with constant scale /i the change of scale is obtained 
by integrating over an insertion of G(x) 

. . . a(a;„))e = jd^x {@{x)0,{x,) . . . 0^{x^)), . (2.3) 

Here 0i , . . . On are local operators and the subscript c at the correlator brackets marks 
connected correlators. 

Assume that a family of renormalizable QFTs is parameterized by renormalized coupling 
constants A*, i = 1, . . . ,N. We assume that an action principle [1] is satisfied. This means 
that for each coupling A* there exists a local operator 0i(x) such that for any set of local 
operators Oi , . . . , On 

-^{Olix,) . . . On{Xn))c = jd^X {Ux)0,{x,) . . . On{Xn))c • (2.4) 

Note that the integrability of the integrand in fl2.3p . fl2.4p assumes the appropriate infrared 
behaviour of the correlators. 

Assume further that the couplings A* can be promoted to local sources A*(x) for the fields 
(f)i{x). The generating functional \nZ then in general depends on the scale factor /i(x) and 
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the sources A*, and the action principle fl2.4l) means that in addition to fl2.ip we have 

(51nZ 



. (2.5) 



A correlation function of the form 

(0,,(X1)0,,(X2) . . .0^„(x2)e(l/i)e(y2) . ■■Qiym))c (2.6) 

evaluated on a flat can be obtained by taking variational derivatives of In Z with respect 
to the sources A* and the metric scale factor /i and then setting the sources and the scale to 
be constant. In a renormalized theory the correlators (12. 6p are distributions. They form a 
basic set of local physical quantities defined in a given QFT. 

In a renormalizable QFT a change of scale can be compensated by changing the couplings 
A* according to (11.11) . By the action principle (12. 4p this implies that Q{x) = /3*0i(x) where 

are the beta functions. This equation should be understood as an operator equation, 
that is, as an equation that holds inside correlation functions (12.61) up to contact terms (i.e. 
up to distributions supported on subsets of measure zero). The use of sources A*(x) and 
non-constant Weyl factor //(x) facilitates bookkeeping of the contact terms. In the presence 
of non-constant A*(a;) and fi{x) one can expand the difference 9(x) — /9*(A(x))(^j(a;) in terms 
of derivatives of the sources and metric [19]. The expansion must by covariant with respect 
to changes of coordinates. This requirement ensures that the contact terms respect the 
conservation of stress-energy tensor. In [19] H. Osborn assumed that this expansion has the 
form 

e(x) - /3V,(x) = ^fi^R2ix)CiX) + d''[W,iX)d^X'] + ^d^X'd^'X^C.iX) (2.7) 



where 



fi^R2ix) = -2d^d''\nfi{x) {2.i 



is the two-dimensional curvature density. Note that in (12. 7p C, Wi and Gij are functions of 
A evaluated on A*(x) that depend on x via A*(x) only. Effectively equation (12.70 gives a local 
version of renormalization group equation. Using the Wess-Zumino consistency conditions for 
the local renormalization group transformations (12.71) H. Osborn derived a gradient formula 

m 

diC + gijP' + bijp' = (2.9) 

where c and Qij are the Zamolodchikov c-function and metric [11] defined in terms of two- 
point functions as 

c = 47r2 (^x^'x^x'^x'^ - x^g'^^x'^x'^ - ^x^x'^g'^'x^^ { T^^(x) r„^(0) )c^^^^^^^ (2.10) 

=67r^A-^(0,(x)0,(O)),y^l^l^^ (2.11) 

where is a fixed arbitrary 2-d distance. The tensor bij is an antisymmetric two-form 
that can be expressed as 

bij = diWj - djWi , Wi = 37i j£xx^9{l- A\x\){(f)i{x)e{0))c (2.12) 

where A is the same mass scale used in the definition of c and g^j. The most restrictive 
assumption in [19] appears to be the form of expansion (12.70 . The fact that the expansion 
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does not go beyond the second order in derivatives suggests a certain power counting prin- 
ciple. Such a principle could be provided in the vicinity of an ultraviolet fixed point by the 
standard power counting arguments for renormalizability. Even with such a counting prin- 
ciple the expansion (12.71) is too restrictive. Thus it omits terms of the form 9^AVf (x) where 
J^{x) are local vector fields which can be prescribed engineering dimension 1. Such terms 
in the scale anomaly can be generated by near marginal perturbations near fixed points. In 
particular they are present in generic current-current perturbations of Wess-Zumino-Witten 
theories |27j. Another class of theories for which (12.71) is too restrictive is general nonlinear 
sigma models. In this case one needs to allow the quantities C, Wi and Gij in (12.71) to have 
a non-trivial operator content. The case of sigma models was covered separately in [12] (see 
also [12], [H], [15], [T7], [T8] and references therein). It was shown that a gradient formula 
analogous to (12. 9p can be derived provided a sigma model integration measure with certain 
properties exists. In the present paper we will go beyond Osborn's UV assumptions allowing 
for an arbitrary local covariant expansion with operator- valued coefficients replacing (12. 7p . 
Making instead assumptions about the infrared behaviour we derive a general formula 

dic + {gij + Agij)(3^ + hj(3^ = . (2.13) 

The metric correction Agij is constructed via two point functions of 0j with the currents 
Jj{x) arising from the expansion generalizing expansion (12.71) (see formulas (18.21) . (I8.15P ). 
Alternatively Agij can be expressed via 3 point functions with the pure-contact field D{x) = 
<d{x) — j3{x) (formula (18.30 ). Formula (12.130 is derived under two separate assumptions on 
the infrared behaviour. The first assumption is that the action principle (12. 4p holds for 
one and two point functions of operators (pi that assumes that these functions are at least 
once differentiable. This ensures in particular that the c-function is once differentiable. The 
second assumption is that for any vector field J^{x) we have 

lim \x\'{J^ix)To.p)c = 0. (2.14) 

This condition is equivalent to requiring that the the long distance limit of the QFT does 
not exhibit spontaneously broken global conformal symmetry. (Recall that at fixed points 
special conformal symmetry requires T{z) to decay at infinity as |^|~*^.) As a simple example 
in section [9^2] demonstrates, this condition is essential. If in a scale invariant theory the global 
conformal symmetry is broken via boundary conditions at infinity the value of the central 
charge may vary with moduli. 

Our considerations include the nonlinear sigma model case. We thus show that in order 
to have a gradient formula we may replace the somewhat obscure technical assumption on 
the measure given in [TB] by a more conceptually clear assumption on the stress-energy 
tensor behaviour (12.141) which we show to be a necessary assumption in section 19.21 A 
question remains, of course, how one can check whether our infrared conditions hold in any 
given theory. Since in the nonlinear sigma model the expectation values of diffeomorphism 
invariant local operators are believed to be free of perturbative infrared divergences they 
must be analytic in the couplings ([21], [25]). This means that the first infrared assumption 
can be controlled in perturbation theory. It is less clear to us whether one can control the 
infrared behaviour of T^,y perturbatively. We are planning to discuss applications of our 
general result (I2.13P to nonlinear sigma models in more details in a separate paper [27j . 
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3 Zamolodchikov's formula 



Zamolodchikov proved in [TT] the following formula 



= -^'g^i^' (3.1) 

where /i is the RG scale, c is the c-function (12.101) and gij is the metric introduced in (12. lip . 
This formula implies that c decreases under the renormalization group flow and is stationary 
exactly at the fixed points, c is normalized so that at fixed points its value coincides with 
the value of the Virasoro central charge. 

Note that the c-function and the metric gij depend on A only through the dimensionless 
ratio A//i, because according to (12.11) and (12. 4p the fields T^y{x) and (pi^x) are densities in 

implying that their 2-point functions take the form 

{(f)i{x) (j)j{0))c = fi'^Fij{fix) , 
(T^,(x)0i(O)), = /i"F^,,,(/ix). (3.2) 

Before we set out to prove the general gradient formula it is instructive to go over a proof 
of formula (13.11) . One way to prove equation (13.11) is to derive alternative formulas for c and 

9ij 

c = - [(fx GA{x){eix)eiO))c (3.3) 

d_ 

dA 

where 

Ga{x) = 3nx^e{l - A\x\) . (3.5) 

These are the formulas for c and Qij that we will use in the proof of the gradient formula. 
Equation (13.10 follows immediately from formulas (13. 3p and (13.40 : 

/^|^ = -A|r = A^r Z^'^; G'A(a;)(e(a;)e(0)), 



g,j = -A^ ld'xGA{x){U^)<PM)c (3.4) 



^^x dA OA 

= -P'g^JP'. (3.6) 

Replacing (0(a;)6(O))c by {P''(f)i{x) j3^(f)j{0))c in the second line is allowed because they 
differ only by a contact term in x, which gives no contribution since the smearing function 
AOGa^x) /dA is supported away from x = 0. 

While formula (13.40 is evidently equivalent to formula (12. lip the equivalence of formulas 
(I2.10p and (13. 3p for c is shown as follows. Combine the special identity in two space-time 
dimensions 

[x^g^^'g''^ - g^'x'^x'^ - x^'x^g"'^ + 2g^"'x''x'^ - x^g^'^'g^f^) { T^,{x)T^p{Q) )e = (3.7) 
with the Ward identity 

9''(r^,(x)T„^(0))e = (3.8) 
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and CPT invariance 

{T^u{x)T^p{0)),= {T^,{-x)T^f,{0))c= {T^p{x)T^,iO))c (3.9) 

to calculate 

d" [(2x'^x"x^-2x2a;V^-a;V"a;'^) ( T^^(a;)T„^(0) ) J = -3x^Q{x)Q{0))c. (3.10) 
It follows from (l3A0D that 

-j £x G'A(x)(e(x)e(o)), 

= Ti j d\9{l - A\x\) 9^ [(2x'^x"x^ - 2x2xV^ - x^g'^'x'^) {T^,{x)T^p{Q) )J 

= 27r' (2a;''a;'^a;"a;^ - x^x^^V^ - x^g'^'^x'^x'^ - x^x'^g'^'x^) ( T^,(x)T„^(0) , ,_ 

/ A|a:|— 1 

(3.11) 

which demonstrates the equivalence of (12.101) and (13.31) . 

4 The proof of the gradient formula (first steps) 

We start by defining a 1-form rj by the equation 

diC + gi,P^ + +ri = (4.1) 

and show that the remainder term r, can be expressed in terms of correlation functions of 
Q{x) and (f)i{x) with the pure-contact field D{x) = Q{x) — /3(x). Infrared behavior of the 
correlation functions will be an important issue, so we introduce an IR cutoff at |x| = L ^ 
and keep track of the error terms. Our assumptions about IR behavior will be designed 
to ensure the vanishing of the IR error in the limit L oo. 
We start out by recasting gijP^ as 

g,,(3^ = Q7r'A-^<p,{x)<P^{0)P^) j ^^^^^^ = Q7r'A-'{<l>,{x)&{0)) y^^^^^^ (4.2) 

which is valid because f3^(f)j{0) differs from 0(0) only by contact terms. This can be further 
rewritten as 

= Jd^xGA{x){<j),{x)e{0)), = y"rf\GA(x)(0,(x)O(O))e (4.3) 

where the scaling property (13. 2p was used on the last step. Finally using (12.31) we obtain 

g,jP^ = fd^y [ d^x Ga{x) ( 0(y) 0,(x) 0(0) . (4.4) 



Formula (14.41) is infrared safe but as we want to impose the IR cutoff systematically, we write 
instead 



g,,P^ + E,= I d'y I d'x Ga{x) ( 0(y) </>,(a:) 0(0) (4.5) 
\y\<L 
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The Ward identity gives the error term 



El = 


[ 


y" J d\ G^{x){T^M 0.(0;) 6(0)), 




J\y\<L 





27r|/V / f/'xGA(x)(T^,(y)0,(x)e(O)), 



/\y\-- 



(4.6) 



which certainly vanishes in the hmit L ^ oo. 

We next turn our attention to the derivative diC. Assuming that c can be differentiated 
with respect to the couphng constants A*, we can write using formula (13.31) for c and the 
action principle (12.41) 



diC 



Sy / (i'x GA(x)(0,(y)e(x)e(O))e. 



Again, we regularize in the IR as 



d'tc = - \ d^y I d^x G'a(x) ( 0,(y) Q{x) 6(0) ), 
y\<L 



(4.7) 



(4i 



Formulas (14.51) and (14. 8 p can be combined to get 



<9f c + g,,[3^ + Ei= d'y d'xG^ix) ( Q{y) 0,(x) 9(0) - Uv) ©(0) )c 

J\y\<L J 

= [ d'y [ d'xGAix) { IPiy) + D{y)] 0,(x) 6(0) - 0,(y) [P{x) + D(x)] 0(0) 

J\y\<L J 

= -6f/^ + / d'y [ d-'xCAix) {Diy)<P,{x)e{0) ~ <P,iy) Dix)eiO)), (4.9) 

J\v\<L J 



'\y\<L 

where we have introduced the 2-form bj". 



y\<L 



d'y / d'xGAix) {(j),{y) 0,(x) 6(0) - 0,(y) 0,(x) 6(0) ), 



Equation (14.91) can be written as 
with 



\y\<L 



d^y / d^x Ga{x) { Uy) D{x) 6(0) - D{y) 0,(x) 6(0) ), 



(4.10) 

(4.11) 
(4.12) 



Equations (14. lip . (14.120 are the main results of this section. We will later show that under 
our assumptions on the infrared behaviour the limits 



diC = lim 9f c , 

L^oo 



bij = lim 



(4.13) 



exist. The error term Ex goes to zero as L ^ oo. The remainder term rf is expressed via 
correlation functions involving the pure-contact field D[x). In order to investigate this term 
we develop a sources and operations formalism for calculating correlation functions of D{x). 
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5 Sources and operations 



In this section we present a general formalism that allows computing correlation functions 
of pure contact fields using functional differential operators acting on functionals of sources 
and metric. The general exposition is somewhat tedious so for the reader's convenience we 
present the most important ingredients necessary to understand the proof of the gradient 
formula in a separate subsection 15.21 



5.1 General formalism 

So far we have introduced the fields 4>i{x) as operators conjugate to the coupling constants 
A* that parameterize a renormalizable 2D QFT. It will be convenient to assume that the 
set (pi is complete in a given class of fields which we denote by JF. The class of fields can 
be a complete set of spin-0 relevant and near marginal fields. We could define such fields 
without a reference to a particular fixed point by requiring that the corresponding coupling 
constant belongs to some family of renormalizable theories with finitely many couplings 
(there are finitely many couplings for which is not identically zero). This will not work 
for the nonlinear sigma models, for which the set of couplings is infinite, but in that case 
we could talk about near-relevant and near-marginal couplings using the engineering scaling 
dimensions introduced via free fields. As yet another possibility we could assume that the 
set {(pi} spans all spin-0 local fields and work with a Wilsonian RG. We will keep the class 
of fields unspecified throughout this section assuming only that JF is closed under RG 
the precise sense of which we will discuss below. In general a field 0{x) is defined via its 
distributional correlation functions with other fields. If 0{x) G J-' the completeness of {(pi} 
means that there are unique coefficients such that the field 0{x) — 0'^(pi{x) has vanishing 
correlation functions with all fields from JF inserted away from x. The field 0{x) — 0'^(pi{x) 
is thus a pure contact field, that is its correlation functions are distributions supported on 
a subset of measure zero in x. We can define ordinary fields 0{x) as fields for which the 
correlations of 0{x) — 0'^(pi{x) are zero as distributions. This means that the distributional 
correlation functions of such fields are obtained from those of the fields (pi{x) by contracting 
them with the appropriate coefficients O*. 

Whatever T we choose it is essential that the trace of stress-energy tensor can be expanded 
in these fields: 6(x) = /?*(A)0j(x). It is worth noting that the set (pi may include total 
derivative fields. Although the correlation functions are independent of the corresponding 
coupling constants the beta functions may be non-trivial and total derivatives may thus 
contribute to 6(x). Let us further introduce sources A*(x) for all fields (pi{x) so that the 
generating functional InZ depends on these sources and the metric scale factor with 
equations fl2.2p and fl2.5l) satisfied. This means that (pi{x) and Q{x) are represented by 
functional derivatives 

5 5 
'^^^") = 5A^' ®^")=^(")m^) ^'-'^ 
which we chose to denote by the same symbols. The action of these functional derivatives on 
InZ generates distributional correlation functions (12. 6p . To facilitate the use of differential 
operators in computing correlation functions we introduce a shorthand notation 

)) = \nZ (5.2) 
(( = restriction of functionals to constant sources and flat 2-d metric (5.3) 
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so 

(x,) ■ ■ ■ Q{y,) ■■■)) = ( 0., (x,) ■ ■ ■ Q{y,) ■■■), (5.4) 

where, on the left hand side, the and Q{x) are functional differential operators (15. ip . 

while on the right hand side they are fields. 

Define operations 0{x) to be first order local differential operators acting on functionals of 
the sources and 2-d metric. The word local here means that the coefficients of the functional 
derivatives in an operation given at x can depend only on the values of X{x), fi{x) and 
finitely many derivatives thereof. An ordinary field 0{x) = O*(A)0i(x) is naturally assigned 
an operation 0{x) = O*(A(a;))0j. Operations of this form we will call ordinary. An arbitrary 
operation 0{x) gives rise to an ordinary field denoted 0{x) via 

( Oix) 0,, (xi) ■ • ■ e(yi) ■■■)c={{ Oix) 0,^ (xi) • ■ ■ e(yi) • ■ ■ )) . (5.5) 

Although the above formula specifies distributional correlation functions containing only a 
single 0_{x) it defines uniquely the coefficients O* in 0_{x) = 0^(j)i{x) and thus in principle 
fixes the correlation functions containing arbitrarily many 0_{x). The ordinary operation 
0^(f)i{x) corresponding to 0_{x) will be denoted by the same symbol 0_{x). Define pure- 
contact operations 0{x) as operations satisfying 0_{x) = 0, i.e., 

{{O{x) = (5.6) 

Then 

((0.,(xi)-- -6(1/1) ■■■0(0))) = (([0,^(xi), 0(0)] ■■■e(i/i) ■■■)) + ■■■ 

+ ((0,(a;i)^^^[e(yi), 0(0)] •■■)) + ■■■ (5.7) 

is a sum of contact terms. 

We would like now to construct an operation for a given operator that can be used in 
computing its correlation functions from InZ. Since we know how to do this for ordinary 
operators it suffices to solve this problem for a pure contact field. Let 0{x) G be a pure 
contact field that does not explicitly depend on A* that is [di,0{x)] = 0. Then we can 
construct a pure contact operation 0{x) by requiring 

(( 0,, (xi) ■ ■ ■ 6(2/0 ■ ■ ■ 0{x) )) = (0,, (xi) ■ ■ ■ e(yi) ■ ■ ■ Oix)), . (5.8) 

This essentially fixes 0{x) because in physical correlators singularities appear only when 
some of the insertions coincide. The only ambiguity in 0{x) is operations annihilating \nZ. 
Any choice however suffices for practical purposes. With this definition given an arbitrary 
operator A{x) G its correlators with the fundamental fields (pi^i^Xk), ^iyi) can be computed 
using the ordinary operation A{x) = A^(f)i and the contact operation 

A{x) = [A^](x) (5.9) 

according to 

(0,,(xi) ■ ■ ■ 0(2/i) ■ ■ ■ A{x)), = (( 0i,(xi) ■ ■ ■ 0(2/i) ■ ■ ■ [A{x) + A{x)] )) . (5.10) 

In the above correlation function the contact terms proportional to 5{x — Xj^) are essentially 
fixed by the action principle (12.41) . The extra contributions arising from the explicit depen- 
dence of the coefficients on A-^ 's are accounted for by commuting the operation A{x) to the 
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left. Similarly the contact terms proportional to S{x — yi^) are fixed by the change of scale 
equation fl2.3p . All contact term contributions proportional to derivatives of delta functions 
are obtained by commuting the pure contact operation A{x) to the left until it annihilates 
((• 

Consider now the operator 6(x). Assuming, as we agreed before, that 0(x) = = [3{x) 
in the operator sense means that = and the field D{x) = 9(x) — P{x) is pure contact. 
As 9(x) does not explicitly depend on A* {fid/dfi and d/dX^ commute) we can define a pure 
contact operation V{x) in accordance with the general rule (15. 9p . (15.101) . The field G(x) is 
special in that it is represented by a variational derivative (15.11) . This implies that 

(( 0,,(xi) ■ ■ ■ e(yi) ■ • ■ [e(x) - p{x) - v{x)] )) = (5.11) 

that can be written more succinctly as a first order functional differential equation on the 
generating functional 

[Q{x) - f3{x) -V{x)]\nZ = 0. (5.12) 

Knowing the pure contact operation V{x) the correlation functions of D{x) with any number 
of 4>i{x) and Q{x) can be calculated as 

( Dix) 0,,(xi) . . . e(yi) . . . )e = (( (e(x) - pix)) 0,,(xi) . . . e(i/i) ■ ■ ■ )) 

= (( mx) - (3{x)), K (^i) • • • 0(yi) ...])) + (( 0n (^i) • • • 0(Z/1) • • • T^{^) )) 

= {{ [<t>^M) ■ ■ ■ Qivi) m^) + m - <s>ix))] )) 

= (( [0,,(xi), Vix)] . . . Qiy,) ...)) + ... + (( 0,^(xi) . . . [0(i/i), P(x)] ...)) + ... 
+d,J'6ix - xi) ( 0,(xi) . . . e(yi) •■■), + ... (5.13) 

where equation (15. lip was used on the second line, {{T>{x) = was used on the third line 
and 

[0,,(xi), P{x)]=6{x-x,)d,,P'(P,{x,). (5.14) 

was used on the last line. 

The form of V{x) is constrained by 2-d covariance and locality. In general it can be 
written as an expansion in derivatives of the sources A* and covariant derivatives of the 
curvature with coefficients being ordinary operations. It is interesting to consider additional 
restrictions on V{x) from power counting rules. We will distinguish two such rules which we 
call a loose power counting and a strict power counting. In both cases the expansion of V{x) 
goes only up to two derivatives in the sources and metric. In the loose power counting rule 
the coefficients can have a nontrivial operator content. Explicitly in this case we can write 

Vix) = ^fi^R2{x)C{x) + d,X'{x)Jt{x) + [W,{x)d,X] + U.Xd^'X^G^.ix) (5.15) 

where C{x), Wi{x), Gij{x) are ordinary spin-0 fields, and (a;) is an ordinary spin-1 field, 
and where the 2-d curvature is given by 

jj'^R^ix) = -29^9^1n/i(x). 

Two comments are in order here. Firstly, notice the appearance of vector fields Jf (x) in 
the expansion. As we defined operations only for spin-zero fields to accommodate fields and 
operations of nontrivial spin we need to introduce new fundamental fields and new sources for 
those fields. While used to obtain distributional correlation functions involving operators of 
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nontrivial spin such sources are always set to zero in the end of a computation. The operation 
V{x) does contain terms proportional to the tensorial sources and their derivatives. However 
our proof avoids using the explicit form of such terms and we will not introduce the tensor 
field sources explicitly not to clutter the computations. Nevertheless the operations like 
Jf (x), when appear, should be understood in this sense. 

Secondly, notice that in the power counting scheme used the operators C(x), Wi{x), 
Gij{x) must have dimension near zero. This means that, using the fixed point language, 
we allow for slightly irrelevant terms to appear in V{x). This is a common consideration 
used for general nonlinear sigma models [5]. The loose power counting thus accommodates 
perturbative nonlinear sigma models. 

If one assumes the UV behaviour is governed by a unitary fixed point, the only dimension 
zero operator is the identity, and the total UV dimension of T>{x) must be strictly 2 then 
the operators C{x), Wi{x), Gij{x) must be all proportional to the identity operator. We call 
this restrictions a strict power counting rule. It applies in a vicinity of a unitary fixed point 
that has a discrete spectrum of conformal dimensions. Under the additional assumption that 
there are no operators (a;) appearing in T>{x) the case of the strict power counting was 
investigated in [19j . 

Finally the case when the only restrictions on T>{x) come from the general covariance and 
locality can be referred to as Wilsonian. We will prove the general gradient formula (12.131) 
in the Wilsonian case. The proof is simplified if we impose loose power counting. We will 
be discussing in parallel how our steps look in that case. 

As a last comment in this section note that due to equation (I5.12p the operation V{x) is 
subject to Wess-Zumino consistency conditions 

[e(x) - (3{x) - V{x), e{y) - P{y) - V{y)] In Z = (5.16) 

which will be exploited in sections 17.3^ 19.31 

5.2 Summary 

Operations 0{x) are local first order differential operators defined on functionals of the 
sources A*(x) and metric. For the fundamental fields 0i(a;) and the trace of stress energy 
tensor Q{x) the corresponding operations are the functional derivatives (15.11) . We intro- 
duced the notation ((Oi(xi) . . . (9(x„))) for a sequence of operations Oi{xi) applied to the 
generating functional \nZ =)) with the result restricted to constant sources and metric (the 
restriction is signified by the symbol ((). 

Given an operation 0{x) one can extract a field from it by restricting it to constant 
sources and metric (15. 5p . The resulting fields are denoted 0{x) and are called ordinary 
fields. Such fields have the form 0_{x) = 0'^(f)i{x). A pure contact operation is an operation 
0{x) for which 0{x) = 0. 

For ordinary fields the distributional correlation functions are completely fixed by those 
of the fields 0j. More generally a given field A{x) equals a linear combination of fundamental 
fields: A{x) = A'^(j)i{x) only up to contact terms. Such contact terms can be stored in a pure 
contact operation A{x) according to (I5.10p . For the trace of stress-energy tensor Q{x) we 
have 6(x) = ^^(pi^x) = P{x) up to contact terms. The corresponding contact terms are 
stored in a pure contact operation V{x). The generating functional satisfies an equation 
[0(x) — l3{x) — T^lx)] InZ = which can be used to compute correlation functions involving 
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the field D[x) = Q{x) —l3{x) according to fl5.13p . The form of V{x) is constrained by locality 
and general covariance. It can be further constrained by a power counting principle. We 
distinguish a strict power counting, which applies to a vicinity of a unitary fixed point with 
discrete spectrum of conformal dimensions, and a loose power counting that is suitable for 
describing renormalizable nonlinear sigma models. For the loose power counting case V{x) 
can be explicitly written as in formula (15.151) . 

6 The Callan-Symanzik equations 

In the operations formalism the Callan-Symanzik equations for correlators involving fields 
(pi{x) and Q{y) can be obtained by integrating equation (I5.13P over x: 



In view of (16. ip the operations V(l)i{x) and T>Q{x) can be interpreted as extra contributions 
to the Callan-Symanzik equations. 
We further notice that 



This follows from the fact that / d^x 0(x) = fid/dfi, J d^y(f)i{x) = d/dX^, and every term in 
T>{y) is proportional to derivatives of A*(|/) and fj.{y). Equations (16.30 imply that there must 
be ordinary spin-1 fields (and ordinary operations respectively) J^(a;) and Jj^{x) such that 




(6.1) 



It is convenient to define the following operations 




(6.2) 




(6.3) 



VQ{x) = -d^J''{x), 



(6.4) 



P0,(x) = -d^Jtix) . 



(6.5) 



If we impose loose power counting, so that V{x) is given by equation ( 15.150 . then 



VQ{x) 
P0,(x) 




so 



J^'ix) 



d^'Cix) 



(6.8) 
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and Ji{x) defined in (16. 5p in general (without any power counting assumptions) coincides 
with the coefficient in the expansion of P(x) based on loose power counting, equation fIS.lSp . 
In general (without any power counting restrictions) since all terms in X>(x) are proportional 
to derivatives of the sources and/or to derivatives of //(x) there exists a scalar operator C{x) 
such that J^'{x) = 9'^C(x). 

The Callan-Symanzik equations (16.11) for the correlation functions at non-coincident points 
(neglecting contact terms) can be now be written as 

(0,,(xi)...e(i/i)...), = /3^— (0,,(xi)-- -6(1/1)... ), + (r0,,(xi)...e(yi)... ), + ... 

+ (0,,(xi)... [-9^ J^(l/i)]... ), + ... (6.9) 

where 

r0,,(xi) = 9,,/?>,(xi) - 9^J^(xi) . (6.10) 

The terms involving the beta functions can be put into the Lie derivative £/3 so that equation 
f l6.9p takes a more succinct form 

d 

[/^^ - ^l3]{<PhiXi) ■ ■ ■6(?/i) 

= { [-d,Jt,{xr)] . . . 6(1/1) ...), + ... + ( 0,^(xi) . . . [-9, J^(i/i)] ...), + ... (6.11) 
7 The proof continued 

We now come back to the proof of the gradient formula which we left at the end of section 
4. We express the remainder term rf of equation fl4.12p in the source-operation formalism. 
The 3-point functions occurring in equation 14.121 can be written as 

( <f>,{y) D{x) 6(0) ), = {{My) ©(0) [6(a;) - /3(a;)] )) + d,[3^6'{y - x){{ 6(0) <f>,{x) )) 

= ((0.(1/) 6(0) Vix) )) + d.P^'iv - x){{ 6(0) 0,(x) )) , (7.1) 
(D(l/)0,(x)6(O)), = ((0,(x)6(O)P(i/))) + 9,/3^52(x-i/)((6(O)0,(i/))), (7.2) 

so 

( 0,(1/) Dix) 6(0) - Diy) 0,(x) 6(0) ),= {{ 0,(i/) 6(0) V{x) - 0,(a;) 6(0) Viy) )) . (7.3) 
Substituting the last relation in equation fl4.12p and using {{V{x) = 0, gives 



d'y Id'xG^ix) {{ Uy) 0(0) V{x) - 0,(x) 6(0) V{y) )) 

y\<L 



d'y Id'xGAix) {{ Uy) [0(0), V{x)] + [0,(1/), V{x)] 6(0) )) 

y\<L 



d'y ld'xGj,{x) ((0,(x) [6(0), V{y)] + [0,(a;), V{y)] 6(0) )) . (7.4) 

y\<L 
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Note that T>{x) is a pure-contact operation, and < A ^ ^ L, so that 

d'y[Q{0),V{y)] = [d'y [Q{0), V{y)] = VQ{0) (7.5) 

\y\<L J 

d^y[Mx)V{y)] = [d^y[Mx)V{y)] = V4>,{0) (7.6) 

\y\<L J 

d'y{{[cP,{y),Vix)] = [d'y{{[<P,iy),Vix)] = [d,,Vix)]=0. (7.7) 

'\y\<L J 

Using these relations in fl7.4p we obtain 



rf = - / rf^yl27r((0,(y)C2(O)))- / rf^x Ga(x) (( 0,(x) ©6(0) )) 
\y\<L 



(7.9) 



d'xGAix){{V(l),{x)Q{0))) {7.i 
where we have defined an operation 

C2{y) = - jd'x\x' [e(y),P(x)]. 
If loose power counting is imposed, P(x) is given by equation fl5.15p . and we have 

C2{y) = - jd'x^x' [-d^d,6\x-y)]C{x) = C{y). (7.10) 
Thus, with loose power counting, 

VQ{x) = -d''d^C2ix) . (7.11) 

We separate rf into two parts 

= r^i+rt2 (7.12) 

d^xGAix) {{V(Pi{x)e{0))) (7.13) 



<2 = -/ d'yl27r{{(j),{y)C2{0)))- d'xGA{x){{<j),{x)VQiQ))) (7.14) 

J\y\<L J 

then investigate each in turn. 

7.1 The IR condition and the sum-rule 

We investigate rj^i first. Our goal is to show that under certain assumptions this quantity is 
proportional to the beta functions. 
We have 

im^ix) = {{V^{x) = (( [-d.Jtix)] (7.15) 

so 

{{V<p,{x) 9(0) )) = -(( d.Jtix) 9(0) )) = -{d.Jtix) 9(0) ), . (7.16) 
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Substituting this expression in equation fl7.13l) we get 

rf^x Ga(x)( 9^ Jf(x) 6(0) )e. (7.17) 



Now we use the technique similar to the one we used in the proof of Zamolodchikov's formula 
(see section [3]). It is straightforward to check that the Ward identity for T^jj{x) implies 

x^d^Jt{x) 6(0) ), = [x^ Jf (x) 6(0) )e - 2x^x^{ Jf (x) T;(0) ), + x^ J^x) T^{0) )J 

(7.18) 

which allows us to perform the integral in equation fl7.17p . obtaining 

rf^, = X, [x'{ J^{x) 6(0) ), - 2x„x'^( J^ix) T^{0) ), + x'{ Jt{x) T^iO) ) J ^^^^^^^ . (7.19) 

What we want however is an expression proportional to Recall that 

G'a(x) = 37rx^^(l - A|x|) (7.20) 



so that 
and 

We write 



Go{x) = Sttx^ (7.21) 
Go{x) - Ga{x) = 3nx^ e{A\x\ - 1) . (7.22) 



rfi = rf^x[G'A(x)-G'o(x)](a^Jf(x)6(0))e 

'|x|<L 



E2 + / d'x [G'a(x) - Go(x)] ( d^J^{x) 0,(0) ),(3^ (7.23) 

'\x\<L 



with 



= / d^xGo(x)(a^Jf(x)6(0))e (7.24) 

J\x\<L 

= 67r2 X, [x^ Jtix) 6(0) )e - 2x,x^( J^x) T^{0) ), + x^ Jf (x) T^{0) ) J /|^.| J7.25) 

We are allowed to replace 6(0) with P^(f)j{0) to obtain equation fl7.23p because Ga{x) — Go{x) 
vanishes for A|x| < 1, so contact terms in the 2-point function make no difference. 

The IR error term E2 will vanish in the limit L ^ cxd if the 2-point functions ( (x) Tq,^(0) )c 
go to zero at large x faster than |x|~^: 

lim |x|3( Jf(x)T,^(0)), = 0. (7.26) 

|a;|— ►oo 

A violation of this IR decay condition would mean that the long distance limit of the quan- 
tum field theory exhibits spontaneously broken global conformal symmetry. Our main IR 
assumption is that such a spontaneous breaking does not take place and equation (17.261) is 
satisfied. 

The condition lim^^-^oo -E2 = is equivalent to the sum rule 



d'xx^d,Jt{x)e{0)), = 0. (7.27) 
Such a sum rule holds for any spin-1 field, given our infrared assumption. 
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7.2 The term rf'^ 



Similarly to (17.231) we want to write rf'2 as an integral over A|a;| > 1 of an expression 
proportional to jS^ . Equation (17. lip , which one obtains when the loose power counting is 
imposed, motivates the following manipulation of equation f l7.14p . Write the first term, using 
equation fl7.2ip for Go{y), 



d'yl27T{{<p,{y)C2{0))) = - / d'y [d,d^Go{y)] {{<Pdy) C2{0))) (7.28) 

y\<L J\y\<L 

then integrate by parts. Equation fl7.14p becomes 



<2 = i?3- / d'yG,{y){{<Pi{y)d,d^C2m) - / d'xG^{x) {{<P,{x)VQ{Q))) (7.29) 

J\y\<L J 

where is an infrared error 

E3 = - [ d'xd,[d''Go{x){{MO)G2{x)))-Go{x)d^{{MO)C2{x)))] . (7.30) 

J\x\<L 

We further rewrite equation (17.290 as 

= Es + E,+ [ d'x [Go{x) - Ga{x)] {{ ^x) VQ{0) )) (7.31) 

J\x\<L 

where 

^4 = - / d\Go{x){{U^)[d^d^G2{0) +VQ{0)])) . (7.32) 

J\x\<L 

The term E^ is identically zero if we assume loose power counting, by equation fl7.1ip . We 
will show in section [73] that in general E4 vanishes as L 00. 

In equation fl7.3ip . the integration variable x is bounded away from 0, so we can substitute 

(( 0,(x) VQ{0) )) = (( I^e(O) 0,(x) )) = ( -d^J^iO) 0,(x) ), (7.33) 

giving 



' 1,2 



E; + E,+ [ d'x [G'a(x) - G'o(x)] ( 0,(a;) 9^ J^(0) . (7.34) 

J\x\<L 

Finally, we will now show that 

d^J'^iO) = P'd.J^iO) (7.35) 
so that r^2 ^Iso becomes proportional to (3\ up to IR errors. 



r,^2 = Es + E,+ / d'x [Ga{x) - Go{x)] { 0,(x) d^J^{0) ),(3' . (7.36) 

J\x\<L 

7.3 The identity 9^J^(x) = pW^J!^{x) 
We want to show that the ordinary field 

K{x) = f3'd^J^{x) - d^r{x) (7.37) 
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is zero, which is to say that all its non-coincident correlation functions vanish: 

{K{x)<p,,{xi)...), = xy^xi,... (7.38) 

In the source/operation formalism, this means that 

{{K{x)(Pi,{xi)...)) = x^xi,... (7.39) 

To show this we first argue that (17.381) is equivalent to showing that 

[V,Dix)])) = }C,{x))) (7.40) 

for some pure-contact operation ICi{x). We then demonstrate that fl7.40l) is a consequence 
of the Wess-Zumino consistency conditions on T>{x). 
It follows from (16. 2p that 

{{K{x) = {{[-d,J^{x)+(3^d,J^{x)] 

= (([pe(x)-/5^P0,(x)] 

= {{[V,e{x)-f3^Ux)] 

= {{P,Di^)] (7.41) 

where D{x) = Q{x) — j3^(f)j{x) is acting here as an operation. This last calculation implicitly 
uses the obvious identity 

{{f3^{Xix))=(3^{X){{ (7.42) 

and its direct implication 

{{[D, (3^{X{x))] = -{{P^{X{x))V = -P^{X){{V = 0. (7.43) 

Now we have 

{{K{x) 0.,(xi) ...)) = {{[V, D{x)] Mxi) ■•■)). (7.44) 
The operation [D, D{x)] commutes with all the 4>i,.{xr) because x ^ Xr, so 

{{Kix) 0,,(xO ■••)) = (( 0.,(xO ...[V, Dix)] )) . (7.45) 

We now need to show that 

((0,,(xi)---[P, Z}(x)])) = x^xi,... (7.46) 

which by (15. 7p equivalent to (17.401) . 

Equation (I7.40p follows from the Wess-Zumino consistency conditions. Recall that we 
have an equation 

0=[D{x) -V{x)])) . (7.47) 
The Wess-Zumino consistency conditions are 

[D{x) - V{x), D{y) - V{y)] )) = . (7.48) 

It follows from 

[e(x),e(i/)] = 0, [e(a:),/3(2/)] = 0, [/3(x),/3(?/)] = (7.49) 
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that 

[D{x),D{y)] = (7.50) 

and therefore (17.471) is equivalent to 

[V{y), D{x)] )) = - {[D{y), V{x)] + [V{x), V{y)]) )) . (7.51) 

The operation ['P(x), 'P(y)] is evidently pure-contact. It also follows from (12. 3p and (17.431) 
that 

[J d'yD{y),V{x)] (7.52) 
is a pure contact operation. Thus integrating equation (I7.5ip with respect to y gives 

[V,Dix)])) = }C,{x))) (7.53) 

where 

/Ci(x) = -[J d'yD{y), V{x)] - [V{x), V] (7.54) 

is pure contact. This completes the proof that at all non-coincident correlation functions of 
f3^dfj_Jj{x) — di^iJ^{x) are identically zero. Therefor^ 

d^J^{x) = /?^a^Jf (x) . (7.55) 



7.4 i?4 is an IR error term 

We owe a proof that the term given by 



^4 = - y d'xG^{x){{<P,{x)[d^d^C2{Q) + VQ{Q)])) (7.56) 

\x\<L 

is an infrared error, that is it vanishes as L — > oo. The argument is a bit tedious, so the 
reader might want to skip this subsection at the first reading. 

We have Note that in general (without the assumption of loose power counting) we have 

[6(0), P(y)] = -d,d^5{y) CM + d,d,d,5{y) G^'^y) + ... (7.57) 

where the omitted terms contain derivatives of delta functions of order 4 and higher. For 
our purposes this expansion can be written more compactly as 

[6(0), D(i/)] = -d,d^5{y)C,{y) + d,d,d,5{y)Cr{y) (7.58) 
where C^^'^iy) is some tensor operation. Formula (I7.58P implies 

5m'9^C'2(0) + P6(0) = d^d,d.,Cl^''\Q) (7.59) 

and therefore 



{(0,(x)[9^9'^C2(O)+P6(O)])) = {d,d,d,C,'''\Q)<P,{x)) + {mx),d,d,d,Crm))- (7.60) 



^It is worth noting that relation 117.551 1 is a generahzation of the Curci-PafFuti relation 1161 known for nonlinear sigma models 
methods similar to those employed in this section one can a 
need this stronger relation in the proof of the gradient formula. 



By methods similar to those employed in this section one can actually prove a stronger relation: Jf'(x) = Jj{x). We do not 



20 



The second term on the right hand side of fl7.60l) vanishes because it is proportional to a one 
point function of a total derivative operator. Thus we obtain 



E4 = -37r j d'xx\<i)i{x)d^d^d^C{ \<d)) (7.61) 

\x\<L 

which exhibits that is a linear combinations of two point functions at separation L. 
Assuming that {(f)i{L)C2, (0)) is integrable at infinity (which is consistent with (C3 (0)) = 
being independent of Aj) all combinations of two point functions entering E4 go to zero as 
L ^ 00. 

8 Conclusion of the proof 

Combining our results for rf^ and rfg, equations fl7.23p and fl7.36p . and substituting in 
equation (17.121) . we get 

rf = E2 + E, + E,+ {Agf^)(3^ (8.1) 

with 

Ag^^ = [ d'x[GA{x)-Go{x)]{(f),{x)d,J^{0)+ (f>,{x)d,Ji'{0)),. (8.2) 

J\x\<L 

The metric correction Agjj can be also written without any direct reference to currents Jf (x) 
using the Callan-Symanzik equations (16. lip 

Ag^^= [ d^x[GA{x)-Go{x)]{Cp-fi^){Ux)^,) (8.3) 

J\x\<L Ofl 

which, using (12. 3p . (12.41) . can be written in terms of integrated three point functions of 
fundamental operators up to IR error terms. 

Equation (14. lip becomes, finally, the IR-regulated gradient formula 

d^c + (g,, + + bf^)P^ + E{L) = (8.4) 

with total error 

E{L) = Ei + E2 + Es + E^. (8.5) 
The L-dependent constituents of the formula are: 

L„ _ / j2„, / ,2 



dfc=- I d'y I d'x G'a(x) ( 0,(y) Q{x) 9(0) ), , (8.6) 

\y\<L 



= / d'y / f/^xGA(x)(0,(y)0,(x) 0(0)- 0,(1/) 0,(x) 6(0) )e, (8.7) 

J\y\<L J 

E, = 2Tiy^y^ j d'x Ga{x) ( U^) ©(0) ) Y,^,^^ , (8.8) 

E, = X, [x'{ Jtix) 9(0) ), - 2x^x^{ J^x) T^{0) ). + x'{ J^x) T^{0) ) J ^^^^^^ , (8.9) 

E3 = - [ d'xd^[df^Goix){{<j)mC2{x)))-Go{x)d^{{<f)mC2{x)))] , (8.10) 

J\x\<L 
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E4 = -3tt J d'xx'{(t)i{x)d^d^d^C£ '{0)) (8.11) 

\x\<L 

and Ag^j is given in (see equations (gSD, dHHD, dMD, dZSS]), dZSHD, fTTBT]) ). 

Now that the infrared regulated formula fl8.4p is derived we can study its L — cxd limit. 
Let us recapitulate our assumptions on the infrared behavior. Firstly, we assume that the 
action principle holds at least for one and two point functions so that the one and two-point 
functions are at least once differentiable. Secondly, the infrared behavior of the stress-energy 
tensor correlators should satisfy (17.261) . The first assumption means that 2,3 and 4-point 
functions involving (f)i{x) or T^y{x) decay faster than when |x| oo. This together with 
formula (I7.26P imply that 

lim E{L) = , 

L— >oo 

lim d^c = diC , 

L— >oo 

lim = (8.12) 

L— >oo 

where bij is given by Osborn 's formulci (nn^ . Note that in showing (KI^ formula ([72SD 
is needed only to argue that E2 vanishes at infinity while the first infrared assumption alone 
suffices to show all other limits. 

Note that although the same set of assumptions implies 

lim Ag^-p^ < 00 (8.13) 

there is no guarantee that the L 00 limit of Ag^j is finite. However, Infrared divergences, 
if present in Ag^^, are orthogonal to the beta function. Therefore they can be subtracted to 
obtain a finite quantity Agij so that the following gradient formula holds 

diC = -{gij + Ag,j + bi,)[5' (8.14) 

where 

Agij = lim [Agf, — subtractions ] . (8.15) 
This completes the derivation of the general gradient formula. 



9 Discussion 

9.1 Contact term ambiguities and scale dependence 

As the proof of the gradient formula uses distributional correlation functions which have con- 
tact term ambiguities one should ask if the formula itself is free from such ambiguities. The 
contact term ambiguities arise from the choice of renormalization scheme and are generated 
by adding to the generating functional finite local counterterms of the form 

\nZ[\g,j\ ^ \nZ[\g,j] + j d'x[f{X)fi'R2{x) + ^c,,{X)d,Xd>'X{x) + . . .] (9.1) 

^The 2-form bij is exact provided Wj defined in 12.121 1 is differentiable. If one relaxes the differentiability assumptions there 
is room for the limit bij = lim^^oo bf^ to exist without Wj being differentiable, in which case bij would be closed but not exact. 
The failure of differentiability of Wj could come from some non-perturbative effects. 
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where /(A) and Cjj(A) are arbitrary function^ (scalar and tensor respectively) and the omit- 
ted terms contain higher order derivatives of the metric and sources. The redefinition (19.11) 
shifts the terms in the renormalization operation V{x). The low order terms shift as 

C{x) ^ C{x)+/3%f{x), (9.2) 
W^ix) ^ Wiix)-d,f{x)-c,jP^{x), (9.3) 
Gij{x) ^ Gij{x) - CpCij{x) . (9.4) 

with all shifts proportional to the identity operator. 

The c-function and the metric tensors Qij, Agij can each be written in a form involving 
two point correlators at non-zero separation only (see formulas fl2.1Up . (12.111) . (18.21) ). Thus 
these quantities are independent of the contact term ambiguities. The 1-form Wi defined in 
(48) changes under (19.11) as 

Wi I— >• Wi — dif (9.5) 

and the antisymmetric form bij thus does not change. Since the redefinition (19. ip is the most 
general on^ the two-form bij is also independent of the contact term ambiguities. 

Another property that we would like to check is whether the quantities we defined depend 
on the scales /i and A only via their ratio fi/A. For the c-function (I2.10p . the metric (12.111) 
and the antisymmetric form (I2.12p this immediately follows from the scaling properties (13. 2p . 
As for the metric correction Agij it may happen that the infrared regulated quantity Ag^j 
contains a logarithmic divergence ~ InL whose subtraction requires introducing a new scale. 
If this happens the subtracted correction will not depend on /x and A via the ratio /i/A only. 
The physical significance of this is unclear to us. 



9.2 The infrared condition: an example 

Here we discuss a simple example that demonstrates the necessity of the infrared condition 
(I7.26P for a gradient formula to hold. Consider a free compact boson X defined on a two- 
dimensional curved surface with metric g^^ by the action functional 

S[R, g,,] = ^ [d'x {X^g^'^d.Xd^X + QX^R,) (9.6) 
87r J 

where A is the coupling constant corresponding to the radius of compactification squared, 
i?2 is the curvature of g^u, and Q is a parameter. Promoting A to a local source A(a;) we can 
define a generating functional 

\nZ[\{x),g^,{x)]= y"[dX] e-^[^(")'^'"'(")] . (9.7) 

For the zero mode integral to be well defined we assume that the theory is defined only on 
a surface with the topology of a plane so that 

j£x^R2 = 0. (9.8) 

and the zero mode integral in ( 19. 7p only yields an overall numerical factor. Note that Q 
cannot be considered as a coupling constant as it does not stand at a local operator. The 

*We assume that these functions are at least once differentiable. 

^The higher order terms omitted in 1 19. Il l do not contribute to the change of Wi. 
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functional integral is Gaussian so the anomaly can be readily computed (e.g. using the heat 
kernel method) with the result 

D{x) = e(x) = ^C{X)^R2{x) + J^{x)d,\ + ^gxxd.Xd'^X + d,{wxd''X) (9.9) 

where 

J'Ax) = ~a-x(x) (9,11) 

The value of wx is essentially scheme dependent. It can be shifted by adding to S" a local 
counterterm J d?x f{\{x))R2{x) dependent on an arbitrary function /(A). In the context of 
nonlinear sigma models such term can be fixed by target space diffeomorphism invariance. 
For the model at hand this gives wx = (SttA)"^. 

We see from (19.101) that while the theory has a vanishing beta function, its c-function: 
c = 127rC(A) has a nontrivial derivative with respect to the modulus A. We can further 
observe that it is the broken global conformal symmetry that is responsible for the breakdown 
of gradient property. The stress-energy tensor on a fiat surface is 

T^u = ^(: d^Xd.X : -% : d-,Xd''X :) + ^{6,^^ - d^d,)X (9.13) 
47r Z 4:71 

It has exactly the same form as the background charge model [26] with imaginary background 
charge. Note that in our theory there is no background charge. Moreover since our theory 
is defined on a topological plane the field X can be taken to be compact with an arbitrary 
radius. The correlation function 

{T{z)Jx,M) = (9-14) 

means that special conformal transformations are broken by the boundary condition at 
infinity, B 

Another way to see the necessity to have a theory defined on a sphere of large radius is in 
the context of nonlinear sigma model. There it is essential for the gradient formula to hold 
(at least in the leading order in the a' expansion) that the zero mode measure includes the 
dilaton contribution corresponding to spherical topology [19j . 



9.3 Bare gradient formula 

Here we will show how the Wess-Zumino consistency condition for the local renormalization 
operation can be used to derive a different gradient formula. The main quantities in the new 
gradient formula are constructed using the anomalous contact terms present in V rather 
than correlation functions at finite separation. For this reason we call it a bare gradient 
formula. As a consequence of that the terms in that formula are defined modulo contact 

® The charge f dz z^T{z) does not vanish at infinity. 
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term ambiguities discussed in section I9.1[ The new formula also suffers from potential 
infrared divergences in the metric. In this section however for the sake of brevity we will 
not introduce an explicit infrared cutoff and our manipulations with integrals will be formal. 
It is straightforward however to introduce such a cutoff with the main result correct up to 
some error terms vanishing when the cutoff is removed. 
Using (17.491) the Wess-Zumino consistency condition 

[{D{X2) - V{X2)), {D{xi) - V{xi))] )) = (9.15) 

can be rewritten a^ 

[Q{x2),V{xi)] - mx^),V{x2)] - [(3{x2)M^i)] 
+p{xi)D{x2) - V{x2)e{x,) + V{xi)D{x2)\ )) = . (9.16) 
Applying to the above equation {{(j)i{y) on the left and integrating over Xi we obtain 

((0,(y)[pe(x2) + {{My) jd''x^p{x^)D{x2))) + {{vMy)D{x2))) 

-l^^{{D{x2)My))c - - X2)d,P'{<P,)) = (9.17) 
where we used the identities 

frfxi[e(xi),P(x2)] =0, {{Uy) jdxi[P'ix2),Vixi)](P,ix2))) =0. (9.18) 

As we know from section 17.31 VQ — (3^'D(f)j is a pure-contact operation. Its field part 
P^d^Jj^ — d^J^ vanishes (is pure contact). Equation (19.171) expresses the contact terms with 
(l)i{y) via the operation T). Integrating the above formula over X2 with the weight {x2 — yY 
and using 

Id'x2 {D{x2)Uy))c{x2 -yf = ^ (9.19) 

we obtaiii^ 
where 



di{C2) = -Hijl3' + CpWi + Qi (9.20) 



Hi, = -G,,-^Jd'yy^[{d,J!;m,iy)) + {d,J.i'mjiy)) 
G., = -7 [d'yy\{mO),VUym, 



4 

4 



W^ = ] ld'yy'{D{y)UO)). 



= \ Jd%y'{d,Jtiy)Qmc. (9.21) 



^Note that this form of the Wess-Zumino condition is linear in 23. This leads to essential simplifications in computations 
and also ensures that terms with tensorial sources in 23 do not contribute to the final gradient formula. 
^Recall that the currents and the metric dj in 115.151 1 are ordinary operations so that {diJ'^) = 0. 
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Note that the tensor Gij is symmetric. This follows from the fact that operations (l)i{y), 
4>j{x) commute. The metric tensor Hij can be also written in terms of integrated correlation 
functions 

1 
4 



d'x {D{x)My)<PM)c - d,p''{My)<Pjmc - djP''{Uy)Mo))c] • (9.22) 



According to our main infrared assumption fl7.26p Qi vanishes and we have a gradient 
formula 

a,c(°) + gifp^ + hfj^(3^ = (9.23) 

where 

c(°) = - W,f5' , g^f = H,, , = diW, - djWi . (9.24) 

The metric H^j potentially suffers from the same infrared divergences as the correction to 
Zamolodchikov's metric defined in (18.151) . We define the finite quantity entering (19.231) by 
subtracting these divergences. 

When loose power counting applies the above quantities can be computed more explicitly 
using (I5.15p . In this case we have 

H^, = G., - ^ ld^y'[{d,Jt{y)<Pm + {d,J^{y)UO))] , (9.25) 

(C2) = (C) and Gij, Wi coincide with the respective quantities defined in formula (I5.15p . In 
the case when the currents Jf are absent formula (19.230 matches with the one obtained by 
Osborn [19]. 



9.4 Dressing transformations 

For any gradient formula 

diC + gi,P^ + bijP^ = (9.26) 

with a symmetric tensor g^j and an antisymmetric tensor bij one can redefine c, hij and gij 

as 

c = c + p'djp^, 

9ij 9ij ^fS^^ij 5 

bij = bij - {di[3c)ij (9.27) 

so that a gradient formula diC = gij(3^ + holds. The tensor Cij above is any tensor on 
the space of couplings that may depend on the couplings and the renormalization scale /i. 
We will refer to redefinitions (I9.27P as dressing transformations. One can show that formula 
(I8.14P is related to formula (I9.23P by means of a dressing transformation specified by 



so that 



4= jd^xG^{x){Ux)<P,{0)), (9.28) 
c = c(°) - p'cfjP^ . (9.29) 
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It is not hard to construct using dressing transformations a class of c-functions that 
monotonically decrease under the RG flow. Such functions c-^ can be defined as 

c^ = -37r jdPxx^f{x^){Q{x)Q{Q))c (9.30) 

where /(x^) is a function such that /(O) = 1, /(x^) decreases fast at infinitjH and 

x^d^f\x^)<Q. (9.31) 
These potential functions satisfy a gradient formula 

d.cf = -{gf+^gUhlW (9.32) 



where 



4 = Jd'xx'[x^d^f{x^)]{Mx)(pjmc (9.33) 



Agi^ = 3n Jd'xxVix') - l]((9,Jf (x)0,(O)) + (9,j;(x)0.(O))) (9.34) 

b{^ = diW^ - djw{ (9.35) 

w{ = 37i jd^xx'^f{x'^){(l)i{x)Q{Q))c (9.36) 
Such smeared c-functions were first considered in 1201. 



9.5 Renormalization group transformation as a flow of couplings 

As one can observe from the form of Callan-Symanzik equations (16. lip the scale transfor- 
mation of correlation functions 

{<Pn{xi)(pi^{x2) . . . e(?/i)e(?/2) • • • )c 

even at finite separation is not fully compensated by the change of couplings A*. In addi- 
tion to changing the couplings according to their beta functions and rotating the fields (pi 
by the anomalous dimension matrices di(3^ the operators 0i(x) and Q{y) each shift by an 
additional total derivative: d^J^{x) and d^J^iii) respectively. If the currents Jf , are not 
conserved these shifts affect the scale transformation of the correlation functions taken at 
finite separation. This signals that more couplings need to be introduced to parameterize 
such additional terms in the Callan-Symanzik equations. Thus to account for the current 
J^iy) it is customary to introduce dilaton couplings that couple to (/)i(x)/i^i?2(x) terms 
in the Lagrangian The generating functional Z depends on these couplings according to 
the functional differential equation 

•^InZ I ^ 5\nZ 

-/i R2[x) . (9.37) 



6X'^{x) 2" '6\^{x, 

The introduction of this new set of couplings is natural if one bears in mind that coupling 
constant redefinitions are responsible for having different RG schemes. To renormalize a 



^An exponential decrease would suffice for all purposes. 
^'^ A completeness of the set ipi is assumed here as discussed in section 5. 
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theory on a curved space one needs counterterms of the form (/)j(x)/i^i?2(a^)- As usual such 
counterterms are defined up to arbitrary finite parts. Changing the dilaton couplings 
accounts for changing the finite parts in such counterterms. (Previously we assumed that 
such counterterms are fixed somehow which amounts to partially fixing the RG scheme. This 
resulted in the extra terms in the Callan-Symanzik equations.) Expanding the operator C{x) 
in flS.lSp as C{x) = P}j(j)i{x) we see that the coefficients P}j can now be naturally interpreted 
as the beta functions for the dilaton couplings. For the loose power counting case the Callan- 
Symanzik equation for correlators of stress-energy tensor takes the form 

d d 
^^-g^{T^,uiyl)Tap{y2) . . . )c = (3'-^^{T^uiyi)Tapiy2) • • • )c + {T^^{yi)Taf3{y2) • • • )c 



+ {T,,{y,)T%{y,) ...), + ••• = {(3^— + (3),—-){TMT^p{y^) . . . ), (9.38) 



where 



= {d^d, - g^,d^d'')C{x) . (9.39) 

We used ( I9.37P and (15.151) to obtain the last equality in fl9.38p . We see that the dilaton 
couplings account for mixings of the stress energy tensor with trivially conserved currents 
r^^(x). With the enlarged set of couplings (A, Ad) the change in scale for correlators of 
stress-energy tensor components (at finite separation) is exactly compensated by the change 
in coupling constants. In particular for the c-function fl2.10p we have 

We can also compute the derivatives of the c-function fl2.10p with respect to the dilaton 
couplings. Using (13. 3p . (19.371) and the identity 

^ ^d'x-^T-, (9.41) 



dX^ J 5\Ux) 



we obtain 

dc d 



dX'r, aAi, 



G.(.)(e(.)e(o)). ^ 2 ji^.GM^)^a^m). • (9.42) 



Integrating by parts in (19.421) . using 

H^^{4>^)=(3^^,{<P,) (9.43) 

and the assumption that 

fd^x (0,(a;)0,(O)) = 9,(0,) = di{(t),) < oo , (9.44) 



we obtain 



- -9?,P' (9.45) 



where 

= 2 jd'x [Go{x) - G^{x)]{<P,{x)d^d^U^)), (9.46) 
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is a symmetric tensor. 

We can further show that the contraction of gradient formula fl8.14p with the beta func- 
tions 13^ gives Zamolodchikov's formula fl3.ip . This boils down to the identity 

3c 

P},^ = P^Ag^.P^ . (9.47) 



(9r 

f^Di^ = 2 Id'x [G^{x) - Go{x)]{Q{x)d,d^C{0)), (9.48) 



Using equations (19.451) . (I9.46P the left hand side of equation (I9.47P can be written as 

_dc 

while for the right hand side we have 

P'Ag^^P^ = 2 Jd'x [Ga(x) - Go{x)]{e{x)P%Jtmc. (9.49) 

The last expression coincides with (I9.48P by virtue of the identity d^d^C{x) = d^J^{x) 
proven in section I7.3[ This identity can be used because the two point function in (I9.49P is 
taken at finite separation. It is not hard to extend the proof of identity (19.470 to a more 
general case not assuming the loose power counting. Formula (19.470 shows in particular that 
the metric correction Agij is necessary to account for the flow of dilaton coupling constants 
when the last ones are present. 

The additional gradient formula (I9.45P together with the main formula (18.141) imply that 
the c-function is stationary with respect to the couplings (A, A/?) at fixed points /5* = 0. The 
inverse follows from the Zamolodchikov's formula (13. ip combined with formula (I9.40p . Thus 
under our main set of assumptions and with loose power counting the stationary points of 
the c-function are in a one to one correspondence with the fixed points. 



10 Final comments 

As we said in the introduction one of the motivations to obtain a general gradient formula 
came from string theory. In regard with potential applications of our result to the problem 
of constructing string effective actions it should be stressed that we worked throughout with 
normalized connected correlation functions while it is the unnormalized and disconnected 
ones which are relevant to string theory. This fact also explains why our results seem to be at 
odds with the conclusion of [15j that the Zamolodchikov c-function does not give a suitable 
string effective action. In the unnormalized correlators the dilaton zero mode 0o contributes 
an overall factor e"^*^" which results in having the same factor in c. Thus stationarity of c 
with respect to 0o implies that c has to vanish at stationary points. This factor and the 
related problem disappear when one builds c out of normalized correlators as we do in this 
paper. 

The aforementioned problem with c prompted various authors to switch to using what 
we call the bare gradient formula which was discussed in section 9.3. The negative side of 
this is that the metric that appear in that formula, being built from contact terms, does not 
have any positivity properties. 

In the present paper we focused on a formal derivation of the new gradient formula and 
discussing its general properties. It would be instructive to illustrate how it works on concrete 
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examples in conformal perturbation theory and nonlinear sigma models. We are planning 
to do this in a separate publication [27]. It is also interesting to understand better the 
implications of the new formula for string theory. We leave this question to future studies. 
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